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We determine the functional behavior near the discrete rotational symmetry axis of discrete 
vortices of the nonlinear Schrodinger equation. We show that these solutions present a central 
phase singularity whose charge is restricted by symmetry arguments. Consequently, we demonstrate 
that the existence of high-charged discrete vortices is related to the presence of other off-axis phase 
singularities, whose positions and charges are also restricted by symmetry arguments. To illustrate 
our theoretical results, we offer two numerical examples of high-charged discrete vortices in photonic 
crystal hbers showing hexagonal discrete rotational invariance. 



PACS numbers: 



I. INTRODUCTION 



Complex scalar solutions of wave equations canpresent 
dislocations similar to those found in crystals The 
essential mathematical property of these complex scalar 
functions in the point or line where a dislocation is local- 
ized is that its phase is increased or decreased in a mul- 
tiple of 27r along a closed curve around it In these 
points or lines, also known as phase singularities, the 
amplitude of the function vanishes and its phase is unde- 
termined. They play an important role in many branches 
of science, such as solid-state physics [3], Bose-Einstein 
condensates (BEC) 13], superfluidity 4], superconductiv- 
ity 0, cosmology [6|, molecular dynamics Q, nonlinear 
optics [1], etc. In the latter case, the study of such sin- 
gularities is often enclosed in a separated branch called 
nonlinear singular optics Q. 

An optical vortex is a complex scalar solution of a 
wave equation defined in a two-dimensional domain char- 
acterised by the presence of a phase singularity [si]. If 
the vortex has a circular symmetry there is only a sin- 
gularity located at the rotational axis. Optical vortices 
with discrete symmetry, or discrete vortices (DV's), have 
been theoretically predicted to exist in inhomogeneous 
periodic media as solitonic solutions of a nonlinear wave 
equation, such as the nonlinear Schrodinger equation 
(NLSE). These periodic media include optically induced 
lattic es photonic crystal fibers or Bessel lat- 

tices [l2|- They have been experimentally observed in 
the former medium [l^, On the other hand, DV's 
have also been predicted to exist as solutions of NLSE 
in self-attractin g B EC in the presence of periodic opti- 
cal lattices [lol.1T5l|. Other solutions with discrete ro- 
tational symmetry and a complicated phase structure 
have been introduced in periodically modulated poten- 
tials^oth in the framework of BEC and nonlinear op- 
tics [H [13, [S H H mi ill . Some of these solutions are 



characterised by the presence of more than one phase sin- 
gularity Finally, other quasi-stationary 
solutions showing discrete symmetry in a homogeneous 
medium, known as necklace beams or soliton clusters 
have been introduced, and some of them show a non- 
trivial phase structure [H, [H [H, [H, [13, HI] . Stationary 
solutions of this kind have been also obtained in an inho- 
mogeneous media such as a photonic lattice [29|, l30|, [Sll ■ 
In this paper, our aim is to study DV's with different 
possible configurations of singularities. We will analit- 
ically obtain the behavior of a DV near the symmetry 
axis to show that they always present a singularity in the 
rotational axis which can be completely characterised by 
discrete group theory arguments. Next, we will show that 
they can present more than one single singularity. The 
positions of these singularities are related according to 
rules arising from discrete group theory arguments. We 
will also provide two numerical examples of DV's with 
more than one singularity to illustrate our results. 



II. THEORY 

To start with, let us introduce first some common def- 
initions which we will be using in the text. If ?A(x) is 
a complex scalar solution of a wave equation defined in 
a two-dimensional domain, x G M^, then the winding 
number 7 of V' along a closed curve F is given by the 
contour integral 7 = 57 ' "^tiere <j) is the phase 

of the complex field V' = IV'le*'^. Let xq be the position 
of a phase singularity of ■0- The topological charge of the 
phase singularity located at Xq is the winding number of 
the complex field ip for the smallest closed curve contain- 
ing Xq. That is, if F^ is a family of closed curves con- 
taining xq parametrized by e such that lim£_>o F^ = xq 
then V = limg^o 57 /p '^4' ' d\. Additionally, one can de- 
fine another quantity, the total angular momentum, as 
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< VI'C^IV'* > / < ^IVj* >, where = {r y. v)^- 

Vortices with circular symmetry can be written as 
i/j{r,9) = g{r)e^^^ , where {r,6) are the polar coordi- 
nates. They present well defined angular momentum 
since TZaip — e^'^V'i where TZa = e^^'" is a continuous 
rotation of angle a £ R and Cz ^ i-^ the generator 
of 0{2) rotation group. Therefore, vortices are eigen- 
functions of the angular-momentum operator, satisfying 
Czip = ^"0- Note that the winding number and the topo- 
logical charge of a singularity are phase-related concepts 
while the angular momentum is a symmetry-related con- 
cept. These circularly symmetric vortices present a sin- 
gle singularity at the origin and fulfill that I = v — 
provided that 7 is calculated in a closed curve that sur- 
rounds the phase singularity. Also, they satisfy that 
g{r) ~ arl'l -|- O (rl'l+^) when r — > Q. 

For a DV, the winding number 7 and the topological 
charge v, i.e. the phase-related concepts, are also well de- 
fined. However angular momentum is not. Nevertheless, 
another symmetry-related concept, the angular pseudo- 
momentum_m, have been defined for discrete symme- 
try media [33]. It has been also demonstrated that this 
quantity is conserved during propagation [s^ . By con- 
struction, the angular pseudo-momentum m completely 
defines the representation of the C„ group to which ip be- 
longs. It has been also shown that m = w = 7 for a DV 
with a single singularity [351]. Nevertheless, to our knowl- 
edge, the relationship among these this three quantities 
has not yet been established for all types of DV's. This 
relationship will permit us to study the number of phase 
singularities associated to a DV. 

Let us show next how the behavior of a DV near the 
symmetry axis can be analitically determined in terms 
of the angular pseudomomentum m. This will permit us 
to establish that there is always a phase singularity of 
charge m located at the symmetry axis. To do this, let 
the function ?/;(r, 9) be a stationary solution of the 2D 
nonlinear eigenvalue equation: 

[La + LNL{\iA)]'4' = f^i^ (1) 

where Lq = —Vf + V{x.) and Li\fL{\t(;\) be invariant 
operators under the C„ point-symmetry group formed 
by discrete rotations of order n around the rotation 
axis. It has been proved in references [13, HJl that 
if ■0 is a self-consistent solution of ([T]) satisfying the 
symmetry condition \ip{r,6 + 2TT/n)\'^ = \ip{r,9)\'^, it 
can be expressed as tpm,a{r,9) = e'^"^^Um.a{r,9), where 
Um,a{''',d) = Um,a{'r,9 + ^), TO is the angular pseudo- 
momentum and a is a band index. The function ^ is 
said to be a symmetric stationary solution of the 2D non- 
linear eigenvalue equation ([T]) if it satisfies the condition 
|V'(r,0 + 27r/n)|2 = |V^(r,0)|2. 

It has been also demontrated that the angular pseudo- 
momentum presents a cut-off related with the order n 
of discrete rotational symmetry of the medium. Par- 
ticularly, it has been shown that \m\ < j for even n 
and \m\ < ^^^^ for odd n [H, [sH]. Under these condi- 
tions it can be proved that the solutions of Eq. ([1]) with 



|m| = 1, . . . , ^ — 1 for even n and |to| = 1, . . . , < 
for odd n cannot be real. It is easy to see that the 
operator L = [Lq + Lnl{\iI!\)] is Hermitian. There- 
fore its eigenvalues are real numbers. Since L satisfies 
L = L* , solutions with m = or to = ^ for even n 
are real (up to a global phase) , since they belong to one- 
dimensional irreducible representations. On the other 
hand, solutions with angular pseudo-momentum differ- 
ent from TO = or to = ^ for even n correspond to com- 
plex solutions belonging to two-dimensional representa- 
tions of the Cn group. Therefore, DV's are characterized 
by |to| = 1, . . . , ^ — 1 for even n and |?7i| = 1, . . . , ^^-^ for 
odd n. 

We leave the study of the solutions with to = or 
TO = ^ for future research. For the rest, i.e., for DV's, 
we will obtain next the mathematical behavior near the 
symmetry axis. Particularly, we will show that if the 
function tp is a symmetric stationary solution of ([1]) and 
satisfies the following mathematical conditions 

1 — ►o 

1. \ip{r,9)\ « \ipo\ + Si}{r,6), where -00 £ R and 
5ip < Vo, 

2. Lnl{\iP\) « Lnl{Ho\) + SLnl, where 
Lnl{\iPo\) e R and SLnl < Lnl{\iPo\) , 
and 

3. V{r, 9) Vo + 5V{r, 9), where Va eR and SV <C 
then 

^p{r, 9) ''S" ar™e™^ + 0(r™+i), 
where a G M. 

From a physical point of view, these mathematical con- 
ditions establish that the amplitude, the nonlinear oper- 
ator, and the potential V have a smooth non-singular 
behavior near the symmetry axis. Therefore, the pre- 
vious conditions are easily satisfied by DV's in discrete 
symmetry media such as the systems mentioned above. 
We will prove next that the mathematical behavior of 
DV's near the symmetry axis depends only on the angu- 
lar pseudo-momentum. 

To demonstrate the above behavior of the func- 
tion ^/j{r,9), we need to express it as ipm.a{r,9) — 
e^"''^Um.a{r,9). The functions Um.a{r,9) satisfy the fol- 
lowing differential equation 

Id TO^ .2to d 19^" 

gj,2 J. Qj, j,2 j,2 QQ J.2 QQ2 ™'" 

+ [V{r,9)+LNLm)] (2) 

Because of the periodic behavior of the wave function 
Um,a{T, 9) and the potential V{r, 9) we can expand them 
in Fourier series in the angular variable as 

k k' 
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and after performing the angular integrals we get the fol- 
lowing set of differential equations for the angular Fourier 
components ^ (r) 



1 d 

r dr 



(m + knY 



E 

k 



(3) 



Assumptions (1), (2) and (3) for the limit r 0, al- 
low us to write V{r,9) « Vq, \tjj{r,9)\ « Voj and 
LNLi\i>\) ~ LnlHtPoD, where Vq, V'o, andLjvid^/'ol) are 
real numbers. Therefore, in this limit, we obtain 



Vk-kir) 



2^ io 

«° S-.kVo + SVk_kir), 



and 



T k — k 



m) 



1 

2^ 



27r 



h.kLNL{\ipQ\) + SL 



k-k 
NL ■ 



Then, the set of differential equations ^ for the an- 
gular Fourier components behave for r ^ as 



dr"^ 

° (Mr; 



Id (m + 
r dr 

J - Vb - Lnl{\'4'o\))u: 



k 



which is the same equation (but in the linear case) 
obtained for circularly symmetric vortices of the form 
4>{r,9) = u{r)e^^^ carrying angular momentum I = 
m + kn. As discussed in the behavior of a solu- 

tion of this type close to the origin is given by u{r) w 
Arl'l + 0(rl'l+^). Therefore, the solution of the previous 
equation behaves in the vortex core region as: 



0{r 



m-\-kn +1 



(5) 



We now analizcd which are the dominant terms in the 
limit r — > in the Fourier expansion of ipm.aiT, S). 



k 



(6) 



In the limit r the dominant term in the expression 
above will be given by that which minimizes the exponent 
of the leading term in the expansion ([51) The behavior of 
Tl{k) — \m + kn\ has to be then analysed for different ni 
and k values in order to find the dominant contribution, 
i.e., that given by the lowest value of ■q{k). 

Let us first consider that, according to the cut-off the- 
orem, |m| < ^ for n even and |m| < ^^^^ for n odd. Be- 
sides, as stated above, for DV's, m ^ and and 



then |m| = 1, ■ • ■ , 2: — l for n even, or |m| = 1, ■ • • , ^^yL 
for odd n. Therefore, it is easy to see that \kn\ > |r7i|, for 
fc| = 1, 2, • • • and, of course, \kn\ < \m\ for fc = 0. With 
this in mind, let us see that 77(0) < ri{k) for all possible 
values of m and k. 

Let us first consider that to > and fc > 0. In this case, 
we have r]{k) = \m + kn\ = m + kn, for fc = 0, 1, . . . , since 
both quantities are positive. Hence, this gives ry(0) < 
ry(fc), for fc = 1, 2, . . . . Next, let us consider that to > 
and fc < 0. For |fc| = 1, 2, • • • , we have 'ri{k) — \m + kn\ — 
\k\n — TO since \kn\ > \m\. Then, this gives r](k) < rj^k + 
1). On the other hand, for fc = 0, we have 77(0) = to. 
Finally, let us show that m < |fc|n — to. For DV's it 
is always satisfied that |to| < ^ in all cases. Then, we 
obtain 2m < \k\n for fc ^ and consequently to < |fc|ri — 
TO. Then, this gives rj{0) < rj{k). 

Now, let us consider that to < and fc > 0. Then, 
we have 7y(0) = |to| and rjik) = \m -\- kn\ = kn — |m|, 
for fe = 1,2, ••• since kn > \m\ and then kn — \m\ is 
always a positive value. Then, this gives r]{k) < ri{k + 1) 
for fc = 1, 2, • • • . On the other hand, since for DV we 
have 2|to| < kn for fc — 1, 2, • • ■ , we obtain \m\ < kn — 
\m\. Therefore, this gives f]{0) < ri{k) for fc = 1, 2, • ■ • . 
Finally, if m < and fc < then we have 77(0) = |m| and 
?7(fc) = |to + kn\ = \kn\ + \m\, for fc = 1, 2, ■ • ■ since both 
quantities are negative. Then, obviously 77(0) < 77(fc), is 
satisfied for |fc| = 1, 2, • • • . 

In conclusion, for all the possible values of m that a 
DV can present, the dominant contribution to the wave 
function arises from the term given by fc = and, con- 
sequently, and according to Eq. ^ and ^ the wave 
function behaves like 

in the limit r — > 0. 

Therefore, it can be easily proved that there exists al- 
ways a phase singularity of charge m located at the point 
Xj. where the rotation axis intersects the 2D plane, i.e. 
the topological charge of this singularity is t; = m. 

Once this has been established, let us go into the rela- 
tionship between m, the winding number, and the topo- 
logical charge of a singularity in depth. The winding 
number of the symmetric stationary solution tp can be 
calculated using any closed curve. Let us consider a 
closed curve F that surrounds the point and let us 
assume that the winding number of the symmetric sta- 
tionary solution ijj is such that 7 ^ to. Then: 

1. There exists always a phase singularity of charge to 
located on- axis, and 

2. There must exist a number of off- axis singularities 
inside the closed curve F fulfilling the condition: 



7 = TO + V 



3 ' 



(7) 



where Vj is the topological charge of the jth phase sin- 
gularity and V is the total number of singularities. 



3. Additionally, there could exist a number of ofF-axis 
phase singularity pairs with opposite charges (vortex- 
antivortex pairs) inside the closed curve T. 

This result is obtained after spliting the path integral 
in the definition of 7, in y + 1 path integrals, all of them 
related with the topological charge of each singularity. 
The integral around the rotational axis will offer topolog- 
ical charge equal to m. The rest must be related with the 
existence of V off-axis singularities or vortex-antivortex 
pairs that contribute with a null net charge to the inte- 
gral. 

Alternatively, let us consider a closed curve F that sur- 
rounds the point for which the winding number of the 
symmetric stationary solution ip is such that j — m. 
Then, using similar arguments, either: 

1. There exists only one phase singularity of charge m 
located on-axis, or, 

2. There exist one phase singularity of charge m lo- 
cated on-axis and a number of off-axis phase singularity 
pairs with opposite charges (vortex-antivortex pairs) in- 
side the closed curve F. 

Hence, we have proved that it can be stated that a 
DV presents a phase singularity of topological charge m 
located at the rotational axis and a number of off-axis 
phase singularities verifying Eq. [71 

Let us show now that the position of these off-axis 
phase singularities of ip is always symmetric with respect 
to the rotation axis. Let us assume that a phase singu- 
larity is located at a point Xg different from the point 
x^, where the symmetry axis intersects the transverse 
plane. Then, the modulus of the function at this point 
vanishes, |V'(xo)| = and the argument is not defined. 
Besides, as stated above, the function ip of a. DV must 
belong to one of the represetations of the group [111, [13] ■ 
Then, if one transforms "0 according to the rotational ele- 
ments of the group, i.e., the transformations C*j — e'*"^, 
t = l,...,n — 1, then the transformed function must be: 



(a) 



(b) 



(8) 



Then, the transformed function and the function itself 
differs only in a fixed phase, i.e. \tpt\ — \ip\ and Z^t — 
+ at, where at = tm'^. 

Let us apply the rotational elements of the group C^, 
t — 1, ... ,n — 1 to the function at the point Xq. The 
transformed function is the function ip evaluated at a set 
of rotated points given by {n^Ot) = {ro,9o + t^), 
t — 1, . . . ,n — 1, where we have used polar coordinates 
and where (rg, ^o) are the polar coordinates of the point 
Xq. According to Eq. (H]), if |V'(xo)| = then {iptl =0 
for t = l,...,n — 1. Consequently, the modulus of the 
function vanishes in this set of points. Let us show that 
there is a phase singularity located in each of these points. 

To do so, let us consider points Xq such that Xq = xq + 



where |x((; | = e with e an infinitesimally small 
) = 9 € [0, 27r]. Taking into account 



quantity and Zx( 

the properties of a phase singularity, the phase of the 
function is increased or decreased in an integer multiple of 
2tt along the circle ^0 obtained after fixing e and varying 
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FIG. 1: Amplitude (a) and phase (b) of a DV with angular 
pseudo- momentum m = ~1 and more than one singularity. 
Red circles in (a) represent the air holes in the photonic crys- 
tal fiber. 



9 from to 27r. Now, let us apply the rotational elements 
of the group C,* = e'*^ , t — 1, . . . , n — 1 to the function 
in all the points in a circle such as ^o- The transformed 
functions are the function tp evaluated at points = 
xt -I- X(j where 9t = 9 + of circles ^t. Taking into 
account Eq. ([5]), the phase of the transformed function at 
each point of the circle and of the function itself at each 
point of the circle ^0 differs only in a fixed number at, 
which is the same for all the points at S,t- Consequently, 
the phase of the function still increases or decreases in 
an integer multiple of 27r along the circles ^t for any e. 

Then, in the one hand, the function vanishes at the 
points Xt, and, in the other hand, the phase of the func- 
tion is increased or dicreased an integer multiple of 27r 
along any circle surrounding these points. This proves 
that, if there exist an off- axis phase singularity, there are 
other n — 1 off- axis phase singularities distributed sym- 
metrically with respect to the rotation axis. 

Then, one can re-writc Eq. ([7]) as: 



K 



J = m + y^^nvk, 

k=l 



(9) 



where Vk is the topological charge of each of the n phase 
singularities related by the symmetry conditions, and K 
is the total number of these rings of singularities. 



III. NUMERICAL RESULTS 

Let us illustrate the previous results with numerical ex- 
amples of optical DV with more than one singularity in 
photonic crystal fibers showing C^v discrete symmetry. 
In Ref. [ll|, [111 I DV's with just one singularity in this 
kind of systems have been introduced. Additionally, the 
cut-off theorem for discrete symmetry media introduced 
in Ref. [s^l states that the only allowed values for the 
angular pseudo-momentum m for a system with discrete 
symmetry of order n = 6, are m = ±1 or m = ±2. In 
Figs. (HI) and ^ we present two stationary solutions of 
Eq. Il]) with self- focusing non-linearity LnlHiPI) = IV'P- 
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FIG. 2: Amplitude (a) and phase (b) of a DV with angular 
pseudo-momentum m = —2 and more than one singularity. 
Red circles in (a) represent the air holes in the photonic crys- 
tal fiber. 

According with the transformations of the Cgu discrete 
symmetry group, these solutions show angular pseudo- 
momentum m — ~1 and to = — 2, respectively. The nu- 
merical calculation of the winding number along a closed 
curve near the boundary of the numerical solutions is 
7 = 5 and 7 = 4, respectively. Then, in accordance with 
Eq. ^ there are V = 6 off-axis singularities in each case, 
and an undetermined number of vortex-antivortex pairs. 
According to Eq. ([9]), there is at least one ring of off-axis 
singularities. A thorough analysis of the phase of both 
stationary solutions shows, first, that there is a singular- 
ity in the symmetry axis. The calculation of the topo- 
logical charge of this singularity is w = — 1 and v — ^2, 
for each case, and therefore equal to the corresponding 



angular pseudo-momentum. On the other hand, there 
are a limited number of points where the off-axis singu- 
larities can be located, i.e., those where the phase seem 
to be undetermined. The calculation of the topological 
charge around these points shows that the singularities 
are located in the points marked with white circles in fig- 
ures ((TJd) and (lib) and present topological charge v = +1. 
Therefore, Eq. (O is fulfilled, since 7 = 5 = -1-1-6 for DV 
with TO = — 1, and 7 = 4= —2 + 6 for DV with to = —2. 
Finally, it is easy to check for both solutions that the po- 
sitions of these off-axis singularities are related according 
to the transformations of the discrete group. 

In conlusion, we have obtained analitically the behav- 
ior of a DV near the symmetry axis. With this result, 
we have been able to establish a general relationship be- 
tween angular pseudo-momentum, winding number and 
topological charge for DV. This rule permits to study DV 
with more than one singularity in any discrete symme- 
try media. Additionally, we have shown that the posi- 
tions of the off-axis singularities are related according to 
symmetry arguments. Finally, the results have been il- 
lustrated with two numerical examples of high-charged 
discrete vortices in a system with Cgu discrete symmetry. 
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